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ABSTRACT 

The significance to which the cosmic microwave background (CMB) observations by 
the satellite COBE can be used to refute a specific observationally based hypothesis 
for the global topology (3-manifold) of the Universe is investigated, by a new method 
of applying the principle of matched circle pairs. 

Moreover, it is shown that this can be done without assuming Gaussian distribu- 
tions for the density perturbation spectrum. 

The Universe is assumed to correspond to a flat Friedmann-Lemaitre model with a 
zero value of the cosmological constant. The 3-manifold is hypothesised to be a 2-torus 
in two directions, with a third axis larger than the horizon diameter. The positions 
and lengths of the axes are determined by the relative positions of the galaxy clusters 
Coma, RX J1347.5-1145 and CL 09104+4109, assumed to be multiple topological 
images of a single, physical cluster. 

If the following two assumptions are valid: (i) that the error estimates in the 
COBE DMR data are accurate estimates of the total random plus systematic error; 
and (ii) that the temperature fluctuations are dominated by the naive Sachs- Wolfe 
effect; then the distribution of the temperature differences between multiply imaged 
pixels is significantly wider than the uncertainty in the differences, and the candidate 
is rejected at the 94% level. 

This result is valid for either the 'subtracted' or 'combined' Analysed Science 
Data Sets, for either 10° or 20° smoothing, and is slightly strengthened if suspected 
contaminated regions from the galactic centre and the Ophiuchus and Orion complexes 
are removed. 

Key words: cosmology: observations — cosmic microwave background — galaxies: 
clusters: general — cosmology: theory 



1 INTRODUCTION 

1.1 Friedmann-Lemaitre models: curvature and 
topology 

At the end of last century, Schwarzschild (1900, 1998) dared 
to suggest 'a few remarks' which require 'a total break with 
the astronomers' deeply entrenched views'. He suggested (i) 
that geometers' abstract concepts of non-Euclidean geome- 
tries could correspond to the real Universe, but also (ii) that 
non-trivial topology could also be required in order to rep- 
resent the space in which we live. Since then, the physical 
theory (general relativity) linking his first suggestion, curva- 
ture, to the stress-energy tensor has encouraged considerable 
observational efforts in measuring (directly or indirectly) the 
curvature of the observable Universe. 

On the other hand, mathematical difficulties and the 
lack of an accepted quantum gravity theory probably explain 



why less observational effort has been directed at following 
up the second of Schwarzschild's two remarks that have 'nei- 
ther any real practical applicability nor any pertinent math- 
ematical meaning', i.e. measuring the topology of the Uni- 
verse. Since the curvature estimates strongly favour a non- 
positive curvature, without measuring the topology of the 
Universe it is not possible to answer fundamental questions, 
such as whether or not the Universe is finite. (Measuring 
topology does not guarantee an answer, but not measuring 
it guarantees not having an answer, except if the curvature 
is positive.) 

Both curvature and topology have been put in a rela- 
tivistic context in the Friedmann-Lemaitre models (de Sitter 
1917; Friedmann 1924; Lemaitre 1958). 

Mathematical developments in three-dimensional geom- 
etry and topology (Thurston 1982, 1997) now abound, and 
theoretical ideas of the generation of topology at the quan- 
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turn epoch are developing rapidly (Madore & Saeger 1997; 
Carlip 1998; Ionicioiu 1998; Dowker & Surya 1998; Dowker 
& Garcia 1998; Rosales 1998; e Costa & Fagundes 1998; 
Gibbons 1998; Ratcliffe & Tschantz 1998). 

Recently, considerable work in observational meth- 
ods, constraints and candidates for the global topol- 
ogy of the Universe has also commenced. These include 
(i) cosmic microwave background (CMB), i.e. essentially 
two-dimensional, methods (see Sect. 1.2) and (ii) three- 
dimensional methods (Fagundes & Wichoski 1987; Fagun- 
des 1996; Lehoucq et al. 1996; Roukema 1996; Roukema & 
Edge 1997; Fagundes 1998; Lehoucq et al. 1998; Roukema 
& Blanlceil 1998). 

Work has also been done considering how simulations 
of large scale structure (the galaxy distribution) are affected 
by the usual assumption (for numerical rather than physical 
reasons) that the Universe is multiply connected (Farrar & 
Melott 1990). 

For theoretical reasons, interest is probably highest in 
the hyperbolic manifolds. However, the observational argu- 
ments tending towards a flat, cosmological constant domi- 
nated universe (Fukugita et al. 1990; Fort et al. 1997; Chiba 
& Yoshii 1997; Pcrlmutter et al. 1999) may bring interest 
back to more thorough analyses of fiat, multiply connected 
Friedmann-Lemaitre models. Moreover, (i) for a cosmologi- 
cal constant to be observable, fine-tuning is necessary if an 
inflationary scenario is accepted, which may also fine-tune 
the topology scales such that r m j < rgLg ; and (ii) the un- 
known physics required to explain Ao > might also explain 
the topology of the Universe. ' 

Hence, the flat models should not be neglected. 

For reviews on cosmological topology, see Lachieze-Rey 
& Luminet (1995); Luminet (1998); Starkman (1998) and 
Luminet & Roukema (1999). 

The 2-D and 3-D methods have different underlying as- 
sumptions on the physics of regions of space or of objects, 
so systematic errors are likely to differ between the two ap- 
proaches. Both have advantages and disadvantages. For a 
list and discussion of both 2-D and 3-D methods, see Ta- 
ble 2 of Luminet & Roukema (1999). A discussion restricted 
to the 2-D methods follows here. 



1.2 Two-dimensional (CMB) methods for testing 
global topology 

The 2-D methods cited above either (i) model (analytically 
or by simulations) the perturbation field in the universal 
covering space''" of compact universes and compare the ran- 



* The shortest distance between two topological images of an as- 
trophysical object is termed the injectivity diameter, 2r m :, and 
can be thought of the 'smallest size' of a small universe. The dis- 
tance to the surface of last scattering (hereafter, SLS) is labelled 

r SLS- 

T The dimensionless cosmological constant is defined Ao = 
c 2 A/(3H^); in the notation of Peebles (1993), A = il\. 
•f The universal covering space can be thought of as an 'appar- 
ent' universe composed of multiple adjacent copies of the physical 
universe. It is a useful, though unphysical, construct for calcula- 
tions. The physical universe itself (one 'copy') is referred to as 
the fundamental polyhedron. 



dom realisations (or statistics) with CMB data (Stevens et 
al. 1993; Sokolov 1993; Starobinsky 1993; Fang 1993; Jing 
& Fang 1994; de Oliveira Costa & Smoot 1995; de Oliveira- 
Costa et al. 1996; Levin et al. 1998a; Bond et al. 1998) or (ii) 
analyse CMB data directly for multiple topological images 
(see Sect. 2.2; Cornish, Spergel & Starkman 1996, 1998b,a; 
Weeks 1998). [The method of Levin et al. (1998b) probably 
falls into the latter category, but is not discussed further.] 



1.2.1 Advantages and disadvantages in the perturbation 
simulation approach 

The fact that the former method (i), (which has been applied 
principally to Euclidean models, i.e. toroidal or hexagonal 
prismical models^) makes statements about an ensemble of 
universe models can be interpreted as either an advantage 
or a disadvantage relative to the latter method (ii). The 
simulational methods (a) test consistency of observational 
and model perturbation statistics, based on (b) statistical 
properties of primordial perturbations, while the identified 
circles principle (Sect. 2.2) directly tests whether or not the 
SLS is consistent with multiple topological imaging, using 
the actual temperature perturbations measured. 

(a) The problem with a result about perturbation statistics 
is that the confidence level with which a particular multiply 
connected model can be rejected is only a statement about 
the rarity of perturbation properties required to be consis- 
tent with the observations. 

On the contrary, the circles principle enables rejection of a 
model by showing that temperatures at supposedly multiply 
topologically imaged points are inconsistent. That is, it can 
show that no set of perturbations would allow the model to 
fit the data. 

(b) In addition, both theoretical and observational justifi- 
cation at large scales of the perturbation statistics assumed 
for the simulational approach are generally made under the 
assumption of simple connectedness. It could be considered 
inconsistent to test multiple connectedness if properties ex- 
pected under the assumption of simple connectedness were 
assumed. 

The theoretical motivations for Gaussian amplitude dis- 
tributions and a P(k) oc k power spectrum are unlikely to 
be valid at scales approaching 2r m j and r+. 

That is, either for a hyperbolic or for a flat, Ao ~ 0.7 
metric to be presently observable (irrespective of topology) , 
inflation needs some degree of fine-tuning [which can partly 
be provided by the ergodicity of geodesies in the former case, 
Cornish, Spergel & Starkman (1996)]. For a flat, Qo — 1 
multiply connected model to be observable, even more fine- 
tuning is needed."' 

So, either flat, observably multiply connected models 
can be rejected for purely theoretical reasons, or else they 
can be tested observationally accepting that perturbation 



3 Bond et al. (1998) applied this same method to two hyperbolic 
models, but using correlation functions instead of eigenmodes. 
Inoue (1999) has demonstrated calculation of eigenmodes of hy- 
perbolic models. 
" f!o is the present value of the density parameter. 
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statistics may differ somewhat from inflationary expecta- 
tions. 

Observationai motivation for Gaussian amplitude dis- 
tributions and a P(k) oc k 1 power spectrum at the 2r m ; and 
r+ scales is equally lacking — if one is testing a multiply con- 
nected universe, since the only observational justification of 
these properties on large scales is that of COBE data anal- 
ysed assuming simple- connectedness. Moreover, although 
several authors find consistency with Gaussian statistics 
(e.g. Kogut et al. 1995, 1996c; Colley et al. 1996), others 
find indications of non-Gaussianity (Ferreira, Magueijo & 
Gorski 1998; Pando, Valls-Gabaud & Fang 1998), which are 
further discussed in Bromley & Tegmark (1999). 

Apart from the problems inherent to the perturbation 
simulation approach, there are some factors which could be 
added to this approach but have not yet been dealt with in 
the work cited above dealing with flat multiply connected 
models: 

(c) Only three of these eight papers discuss T 2 models, and 
only one compares such models with the COBE observations 
in any detail. 

(d) The integrated Sachs- Wolfe effect (ISW) is not consid- 
ered in any of the flat model papers. [Note that the ISW is 
discussed in Bond et al. (1998) for the hyperbolic case.] 

(e) The possible contribution of redshifts in the SLS radi- 
ation due to peculiar velocities at the SLS (hereafter, 'the 
Doppler effect') is not considered in any of these papers. 
This is expected to contribute less than the naive Sachs- 
Wolfe effect to temperature fluctuations ST/T on COBE 
scales, under the assumption of simple-connectedness. How- 
ever, if this assumption is dropped and r m -. < rgLg is being 
tested, it is difficult to be certain about properties of the 
velocity field on scales near r m j based on present quantum 
cosmology theory. 

A conservative observational approach would be to consider 
statistics which take into account this locally anisotropic ra- 
diation. 

(f) There is also a more local observational problem with 
the work done so far for flat and hyperbolic models. What 
is important for falsifying a small universe hypothesis is the 
significance of individual fluctuations, and whether or not 
these are multiply imaged. 

The significance of individual fluctuations has been studied 
by Cayon & Smoot (1995). This paper suggests that some of 
the most significant hot spots {above the 20 degree galactic 
cut) are galactic contaminants, from the Orion and Ophi- 
uchus complexes. Removal of such non-axisymmetric con- 
taminants has not yet been discussed in the small universe 
literature. 



1.2.2 Advantages and disadvantages in the identified 
circles principle approach 

Although tests based on the identified circles principle 
(Sect. 2.2) avoid the problems (a) and (b) above, they share 
a weakness with the perturbation simulation approach by 
being adversely affected by noise, in particular from detec- 
tor noise and from residual foreground contamination, which 
can make genuinely matched circles look different. 

The ISW and Doppler effects will also provide some 
contributions to CMB maps. The former may be treated 



as foreground contamination. The latter will provide locally 
anisotropic emission, which should be taken into account as 
a modification of Cornish, Spergel & Starkman's identified 
circles principle (Sect. 2.2.2). 

The former is a weakness of applications of the iden- 
tified circles principle, shared by the simulation approach, 
for non-flat models and models with a non-zero cosmological 
constant. The ISW, due to time varying of the gravitational 
potential, is likely to be strong if the Universe is hyperbolic, 
particularly at low redshifts and at the large angular scales 
present in the COBE maps (Inoue 1999). In that case, tem- 
perature fluctuations at low redshifts are mixed in with the 
signal from the SLS, so that an identified circles analysis or 
a simulation analysis restricted to the SLS cannot account 
for the full signal at COBE scales. 



1.2.3 Optimal comparison of multiply topologically imaged 
points 

In view of the above comments, it is clear that for the test- 
ing of flat, small universe models, there are several advan- 
tages if the identified circles principle (Cornish, Spergel & 
Starkman 1996, 1998b, a; Weeks 1998) can be used, since 
it avoids (a) making statistical statements about perturba- 
tions and (b) requiring assumptions about the statistics of 

the perturbations, 'I even though it shares the problem of the 
perturbation simulation approach in its dependence on the 
statistics of detector noise and residual foregrounds. 

There already exists a specific candidate for a small, 
flat universe which has been suggested by observations 
(Roukema & Edge 1997; Roukema & Blanlceil 1998). Be- 
cause the candidate is fixed in astronomical coordinates, this 
offers a good example on which to see if COBE data can 
be used to reject a small universe hypothesis by using only 
the COBE data itself, without making theoretical hypotheses 
about the perturbation properties. 

1.3 Observational rejection of Roukema & Edge's 
(1997) T 2 candidate 

The purpose of this paper is to see if this hypothesis, con- 
sidered as a null hypothesis, can be rejected purely on the 
basis of the COBE observations, of about ~ 300 effectively 
independent pixels at which temperature fluctuations have 
been measured, without making assumptions regarding 
the statistics of primordial density perturbations. 



1.3.1 TheT 2 candidate 

The candidate T 2 3-manifold (Roukema & Edge 1997; 
Roukema & Blanlceil 1998), according to which two galaxy 
clusters at redshifts z ~ 0.4 would be topological images of 
the Coma cluster, is presented in Sect. 2.1. 



I The method of Levin et al. (1998b) can probably also avoid (a.) 
and (b). 

** Since the DMR (Differential Microwave Radiometer) is a dif- 
ferential instrument, the beam-width of 7° per antenna implies 
a width of 10° smoothing in the difference maps; a 20° cut for 
the galactic plane is adopted here; hence about 300 independent 
pixels. 



© 0000 RAS, MNRAS 000, 000-000 



B. F. Roukema 



Apart from the striking fact that these three cluster im- 
ages form a near right angle (in 3-D) and that right-angled 
'toroidal' multiply connected universes are those which have 
been most extensively compared to observations, the rela- 
tive simplicity of the hypothesis makes it ideal for testing 
if observational methods of refutation can really be made 
statistically watertight, in preparation for future candidate 
3-manifolds, which are unlikely to be this simple. 

Methods of refuting the supposed identity of the three 
clusters by their individual properties (positions and veloc- 
ities of galaxies, overall masses, X-ray fluxes), by implied 
multiple imaging of large scale-structure, by predicted posi- 
tions and redshifts of further topological images of the clus- 
ter and by proposals of observations relating to sub-cluster 
structures are discussed at length in Sect. 4.2, Sect. 4.3 and 
Sect. 4.4 of Roukema & Bajtlik (1999). 



1.3.2 How to apply the identified circles principle to reject 
the T 2 candidate 

The basic principle of any method used to constrain or de- 
tect candidates for the 3-manifold is that if the hypothesis 
of trivial topology is not adopted, then in some cases, pho- 
tons can travel between a single, physical object or region 
of space to the observer by several different geodesies of (in 
general) widely differing lengths. That is, the objects or re- 
gions of space are observed at widely differing positions and 
redshifts. 

In the case of the CMB, since the redshift is nearly equal 
for all points on the SLS by definition, only special subsets 
of angular positions on the SLS would be multiply imaged. 

Cornish, Spergel & Starkman (1996, 1998b) realised 
that these special subsets form a well defined and easily 
described set in the covering space: the subsets are pairs 
of matched circles. The explanation for this is subtle, but 
elegant, and is briefly described in Sect. 2.2. 

If the radiation from the SLS is locally isotropic, then 
the observed temperature fluctuations observed around two 
circles of such a pair should be identical, to within the ob- 
servational uncertainty. 

If the values of corresponding points on the circles are 
inconsistent, then the hypothesis can be rejected with a con- 
fidence level based on the estimated uncertainty per resolu- 
tion element, or pixel. 

The assumption of local isotropy seems reasonable for 
the COBE data, due to the large smoothing length (effec- 
tively 10°). At this scale, the effective temperature corre- 
lations should be determined uniquely by the naive Sachs- 
Wolfe effect (Sachs & Wolfe 1967), i.e. by gravitational red- 
shifts depending on the depths or heights of the potential 
wells. 

As long as the gravitational potentials are spherically 
symmetric, i.e. if the equivalents of filaments and great walls 
on a linear scale do not exist at z ~ 1100, then this should be 
a fair assumption. If filaments and walls did exist, then since 
temperatures in two different COBE 'pixels' on correspond- 
ing circles represent averages over nearly two-dimensional 
spherical sections at different angles in the fundamental poly- 
hedron, some modelling of the three-dimensional two-point 
correlation function would be needed in order to calculate 
the probability of these averages being identical. 

For this paper, the averaging over the different spherical 



sections, each of diameter ~ lOOOh x Mpc (for Q. = 1) and 
10° smoothing, is assumed to result in identical values. 



1.3.3 How to treat anisotropic (Doppler) radiation when 
applying the identified circles principle 

The amplitude of the Doppler effect at the SLS, which is 
locally anisotropic, is assumed to be negligible relative to 
the naive Sachs- Wolfe effect. Since the first Doppler 'peak' 
is on scales smaller than those resolved by the COBE DMR 
experiment, the fractional contribution of a Doppler compo- 
nent to an identified circles test of COBE data is probably 
small. 

However, if 'linear' filaments and walls existed at the 
recombination epoch, or if velocity fields related to the r; n j 
scale existed, then this might not be correct. 

It should also be remembered that the intrinsic spatial 
separations in a multiply connected universe can be much 
smaller in the fundamental polyhedron than in the cover- 
ing space, so coherent velocity fields on apparently large, 
but inherently small, scales may generate a Doppler effect 
which would be unexpected under the simple connectedness 
assumption. 

Although the naive Sachs- Wolfe effect is assumed here 
to be dominant and to average out to locally isotropic val- 
ues, an extension of Cornish, Spergel & Starkman's (1998b) 
statistic which should be sensitive to Doppler dominated 
temperature fluctuations is suggested here, just in case the 
Doppler component is more important than expected. 

If the Doppler effect dominates, then for pairs of mul- 
tiply imaged pixels which also have nearly parallel lines-of- 
sight, the values should be nearly equal in absolute value, 
with equal signs if the angle between them is close to 0° 
and opposite signs if the angle is close to 180°. Pixels for 
which the lines of sight are roughly orthogonal provide no 
theory-free constraint on topology if the Doppler effect is 
dominant. 

Hence, a correlation function weighted by the dot prod- 
uct of the two line-of-sight vectors of each pixel pair should 
measure the degree to which the sets of points represent 
identical velocity induced temperature fluctuations. A dot 
product weighted correlation statistic is presented [eq. (4)] 
and demonstrated below. 



1.4 Structure of paper and conventions 

In Sect. 2, the hypothesised 3-manifold is briefly presented 
(Sect. 2.1), the matched circles principle is explained and 
the statistics used to estimate the consistency of multiply 
imaged 'pixels' on the CMB are defined (Sect. 2.2), and the 
COBE data are briefly introduced (Sect. 2.3). Results are 
presented in Sect. 3 and discussed in Sect. 4. Conclusions 
are in Sect. 5. 

The Hubble constant is parametrised here as h = 
.ffo/lOOkm s _1 Mpc -1 . Comoving coordinates are used (i.e. 
'proper distances', Weinberg 1972, equivalent to 'conformal 
time' if c = 1). The metric assumed is flat, i.e., Qq + \q = 1, 
since the candidate 3-manifold is T 2 x R. The case Ao = 
is the main case discussed here. 
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Table 1. The positions of the galaxy clusters Coma, 
CL 09104+4109 and RX J1347.5-1145 shown as J2000 right as- 
cension (a) and declination (<5) in decimal degrees and redshift z. 
[Note: the label 'J2000' in Fig. 1 of Fabian & Crawford (1995) is 
interpreted as a typographical error for 'B1950'.] 



„ SLS of (3,1) 





a 


<5 


z 


Coma 

CL 09104+4109 

RX J1347.5-1145 


12.9969 
9.2289 
13.7918 


27.9807 
40.9428 
-11.7533 


0.0239 
0.442 
0.451 


2 METHOD 

2.1 The T 2 candidat 


e 







Roukema & Edge (1997) pointed out that massive clusters of 
galaxies, dominated by their hot, X-ray emitting gas, should 
be relatively good 'standard candles' for observational con- 
straints on the topology of the Universe. This was, in a sense, 
an update of the 'classical' searches for multiple topological 
images of single objects (Gott 1980). 

Although not a goal of that work, it was noticed that 
among a small number (seven) of very bright clusters used 
to illustrate the principle, three of these form nearly a right 
angle (to within ~ 2—3% depending on Ao) with nearly equal 
(to 1%) arm lengths. This happens to be just what would be 
expected if space is T 3 with two side lengths equal and the 
third unknown, or if it is T 2 x R. An a posteriori estimate 
of the probability of this occurring by chance is difficult to 
make in an unbiased way, so observational refutation is more 
prudent than theoretical refutation. 

Moreover, this 3-manifold is very simple, relative to 
other fiat 3-manifolds or to the hyperbolic 3-manifolds 
[which are considered more likely for several theoretical rea- 
sons, e.g. Cornish, Spergel & Starkman (1998a); e Costa 
& Fagundes (1998)]. The implied positions and redshifts 
of multiple topological images of objects are therefore rela- 
tively easy to calculate and observationally refutable predic- 
tions can be made. A list of the various observational argu- 
ments for and against the candidate, apart from the present 
study, is provided in the discussion section of Roukema & 
Bajtlik (1999). 

The fundamental polyhedron (hereafter, FP) considered 
is therefore a nearly square, long prism. Two axes are defined 
as, L x e x , the vector from the Coma cluster to the cluster 
CL 09104+4109 and L y e y , the vector from the Coma cluster 
to the cluster RX J1347.5-1145. The third axis L z e z is de- 
fined as the cross product of these two vectors and assumed 
to be larger than the horizon diameter (L z > 2Rh). (The 
vectors e x , e y and e z are unit vectors.) The positions and 
redshifts adopted for the clusters are shown in Table 1. 

The values of L x w L y vary slightly from ~ 960h _1 Mpc 
to « 1230h _1 Mpc for < A < 0-9. L z is larger than the 
horizon diameter. 



2.2 The Matched Circles Principle 

The matched circles principle (Cornish, Spergel & Stark- 
man 1998b) is illustrated in Fig. 1 and explained in this 
figure's caption. Topological images of the observer in the 
lower half of the image are redundant (they give the same 
set of identified circles as the upper half). Depending on 




Figure 1. Geometry in covering space showing how an observer 
sees multiple copies of circles on the CMB. Multiple copies of the 
observer which need to be considered are shown as black dots, and 
can labelled (i,j) in units of the size of the Universe in the two 
small, fundamental directions. For clarity, the size of the Universe 
relative to the SLS (surface of last scattering) is larger than that of 
the candidate under consideration in this paper. The copy of the 
observer at (0, 1) is omitted for clarity. An image of the observer, 
for example, at (3, 1), is separated from the observer by a distance 
A. The intersection of the SLS spheres of the two observers is a 
circle, of which two points lie in the plane of the figure. One such 
point is at A. Because the two observers are equivalent, the circle 
as seen by the observer at (3, 1) is equivalent to a circle seen in 
the same direction by the observer at the origin. That is, observer 
(3, 1) looking towards A is equivalent to observer (0, 0) looking 
towards A', and similarly for B and B'. These two multiply imaged 
points in the plane are part of the the multiply imaged circle in 
three dimensions. The radius of the identified circles is rtj. No 
identified circle is a great circle [except if (0, 0) were counted as a 
multiple image of the observer, which would not be useful]. The 



SLS has radius r. 



SLS' 



The dashed circle of radius 2rgTC is the 



maximum extent to which multiple images of the observer need 
to be considered. 



the value of the cosmological constant, Ao, the ratio of rgLg 
to L x ~ L y for the T 2 candidate under consideration is 
6.0 <r SLS /L x < 13.2 for < A < 0-9. This implies about 
three to six times as many circle pairs as shown in Fig. 1. 
Since some circles fall within or mostly within 20° of the 
galactic plane, the number which can be usefully compared 
is smaller than this by about a third. 



2.2.1 Locally Isotropic Radiation: the Sachs- Wolfe effect 

To see if the measured temperature fluctuations are consis- 
tent with the hypothesis of physical identity of the circles, 
the most straightforward approach is to test the hypothe- 
sis that temperatures on corresponding pixels are equal to 
within observational error. 

This hypothesis is tested here by considering the differ- 
ence in the temperature fluctuations in two corresponding 
pixels on matched circles to be a random realisation of a 
Gaussian distribution centred on zero with a width deter- 
mined by the uncertainties of the measurements in the two 
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pixels. By normalising the distribution for each pair to have 
a standard deviation of unity, the full set of pairs of multi- 
ply imaged pixels can be combined to form a large sample 
of a single distribution of mean zero and standard deviation 
unity. 

This standard deviation is a, defined 



a 2 (r) 



(f),-(f)„ 



[A(f)]; + [A(ix)]; 



(1) 



where (-^-) . and (^f-J . are the temperature fluctuations in 
the two (hypothetically) multiple images of a single region 
of space-time. 

The assumption that the mean difference is zero is 
likely to be a good approximation independently of whether 
or not the hypothesis is correct, since the dipole compo- 
nents have been subtracted off the data and little correla- 
tion should exist on such large scales if the assumption of 
simple-connectedness is correct. Nevertheless, it is useful to 
check the mean difference, which is defined 



(?).-(?) 



{[Mf)L 2 +[Mf)];} 



1/2 



(2) 



It is also useful to consider the statistic of Cornish, 
Spergel & Starkman (1998b), which is essentially a two- 
point autocorrelation function normalised by the variance 
per pixel [eq. (2) of Cornish, Spergel & Starkman], though 
the interpretation is slightly less straightforward. This statis- 
tic is 



2(f),(f) : 



(3) 



2.2.2 Doppler Effect 

While it is possible that with the resolution of COBE, the 
Doppler effect averages out from three-dimensional space to 
the one-dimensional corresponding circles in a way which 
makes it negligible, the following variation on eq. (3) is pro- 
posed. This modified statistic is 



Sb. 



opp 



2 (¥)«("),*■*. 



(f);+(f 



(4) 



In a case where two corresponding pixels have orthogo- 
nal line-of-sight vectors, the two temperature values are un- 
related (except if the three-dimensional centre of the pixel 
happens to be at the centre of a spherically symmetric grav- 
itational potential), so a noise value would be contributed 
to S [eq. 3], but nothing would be contributed to Srj opp . 

In the case of corresponding pixels whose lines-of-sight 
are parallel and equal, contributions would be made both to 
S and to Srjopp. In the case of corresponding pixels whose 
lines-of-sight are parallel but opposite, a negative value (or 
anti-correlation) would be added to S, but a positive value 
would be added to Sj^opp- 

For these reasons, if the Doppler effect (or another di- 
rectional effect) were present (even partially) , then the value 



Table 2. Weights for the three frequencies of DMR measurements 
used for the 'subtracted' and 'combined' maps (Bennett et al. 
1992, 1994). 



31GHz 53GHz 90GHz 



subtracted 
combined 



-0.34 
-0.49 



0.82 
1.42 



0.70 
0.18 



of S would be low and noisy in the presence of genuinely 
matched circles, to the extent to which the effect is present. 
On the other hand, Spopp would have a high value, although 
less than the ideal value of unity as would be expected for S 
in the case of locally isotropic emission and a zero intrinsic 
(non-topological) auto-correlation function. 

Although Srjopp is calculated here out of interest, it is 
most likely to be useful in application to future measure- 
ments by the Planck and MAP satellites. 



2.3 OBSERVATIONS 

The COBE DMR observations and results are discussed by, 
for example, Bennett et al. (1994), in which a spherical har- 
monical analysis of the two year results is presented. 

In this paper, the results of four years' data are used, 
made available as 'DMR Analysed Science Data Sets' (here- 
after, ASDS) by the COBE team'T, corrected for galactic 
emission either by the 'subtraction' or 'combination' tech- 
niques of removing synchrotron, dust and free-free emission. 
(Hereafter, the 'subtracted' and the 'combined' maps respec- 
tively.) The subtracted maps have the higher signal-to-noise 
ratio of the two. The weights for the three frequencies used 
to obtain the two maps (essentially those of Bennett et al. 
1994) are listed in Table 2. 

Because of galactic contamination, data between galac- 
tic latitudes of —20° and +20° are not considered. 

Analysis of the significance of individual fluctuations 
(as opposed to their statistical significance) was carried out 
by Cayon & Smoot (1995). These authors noted (Sect. 5 
of Cayon & Smoot) that some regions outside of the above 
galactic cut could also be galactic contaminants. The possi- 
bility that the second, third and seventh hot spots listed in 
Table 2B of Cayon & Smoot could be galactic contaminants 
due to the Ophiuchus complex and the Orion complex is 
considered here. 

Since the data set is oversampled, a smoothing by a 
Gaussian of minimum full width half maximum (FWHM) 
of 10° is necessary to extract a significant signal. In fact, as 
is shown below, the result is stronger if smoothing with a 
FWHM of 20° is considered acceptable. 

More recent analyses of galactic contamination of the 
DMR data have been discussed extensively by Kogut et 
al. (1996a, b). Based on these discussions, alternative ways 
of combining the data from the six DMR antennae in or- 
der to further minimise galactic contamination relative to 
the ASDS maps provided by the COBE team, e.g. using a 
weighted average of the 53 and 90 GHz maps plus a 'custom' 
Galaxy cut, would be recommended for further applications 
of the circles principle using COBE data. 

TT At http://www.gsfc.nasa.gov/astro/cobe/cobe_home.htrnl. 
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Table 3. The values of the statistics d, cr, S and Sj^ or> „, as defined 
in eqs (2), (1), (3) and (4), for the COBE data. 



,S' 



% 



£2E_ 



'Subtracted' Analysed Science Data Set: 

10° FWHM smoothing 

-0.04 2.06 0.00 -0.01 

20° FWHM smoothing 

-0.20 2.71 -0.02 -0.01 

'Combined' Analysed Science Data Set: 

10° FWHM smoothing 
-0.01 1.54 0.02 -0.01 

20° FWHM smoothing 
-0.14 1.93 0.02 -0.01 

20° FWHM; excluding G.C., Orion, Oph loop 
-0.02 2.05 -0.03 -0.03 



3 RESULTS 

3.1 Difference parameters: a and d 

The values of the parameters representing pixel pair differ- 
ences and correlations are listed in Table 3. Depending on 
whether the 'subtracted' ASDS or the 'combined' ASDS is 
used, and on whether smoothing is at the 10° or the 20° 
scale, the distributions of temperature differences in cor- 
responding pixel pairs are clearly wider than the a = 1 
distribution expected if the temperatures were intrinsically 
identical and the differences were only due to random mea- 
surement error. 

Formally, a Kolmogorov-Smirnov test on the list of the 
normalised temperature differences used in calculating a and 
d rejects the possibility that the distributions of these differ- 
ences are drawn from a Gaussian distribution of mean zero 
and standard deviation unity by at least 94%. (The weakest 
rejection in Table 3 is for the 'combined' data set with 10° 
smoothing; the other rejections are stronger.) This is the 
case whether or not the N ~ 1000 pixel pairs are considered 
to be independent, or whether only ~ 150 of these pairs (half 
the number of independent pixels) are considered to be in- 
dependent, fn either case, the large number of realisations 
of the difference distribution results in a precise value of the 
dispersion. 

This result is only valid if the COBE team's analysis 
of the random errors is correct, and if contributions from 
systematic errors are negligible. For example, if systematic 
errors or underestimates of the random errors caused the to- 
tal uncertainties per pixel to be twice as high as the random 
uncertainties presented in the data set, then only the 20° 
smoothed 'subtracted' data set would provide a significant 
contradiction. 



3.2 Matched circles 

What are the actual values of the temperature fluctuations 
along the identified circles? Some examples of these are plot- 
ted in Fig. 2 and the positions of two large pairs of circles 
are shown in Fig. 3. 

From this figure the reader can judge what typical frac- 
tions of the sky are available for comparison of matching 
circles, given the galactic cut of ±20° latitude and whether 
any systematic effects might be present. 



Some of the pairs of circles look impressively similar. 

Since the data is smoothed to 20° ~ 2000h~ 1 Mpc, in 
the case of pairs of circles for which A < 2000h~ 1 Mpc, it 
is clearly possible that some identity between corresponding 
circles can be induced just by the smoothing. For example, 
the temperature fluctuations around the two circles in the 
panel for (i,j) = (0, 2) look very similar, but this is explain- 
able by smoothing as much as by physical identity. 

In the panels for (-6,3), (-5,2), (-4,1) and (-4,2) 
there are peaks which appear to coincide, though the over- 
all slopes vary. In fact, these appear as much stronger, co- 
inciding peaks in the 'subtracted' map. This suggests an 
explanation as a galactic feature. In other panels, such as 
(—6,4), (—6,5) (—6,6) and (—5,5) there are strong peaks 
which are highly in conflict between the two circles. 

In the latter three cases, it could be argued that the 
temperature differences cannot be used to refute the T 2 can- 
didate under consideration, because there could be contam- 
ination from the galactic centre, which may not have been 
subtracted perfectly, or may simply have left residual noise 
even if well subtracted. A large part of the central conflicting 
peaks in these three panels is indeed removed if an exclusion 
zone of 60° around the galactic centre is applied. 

Moreover, one could also remove areas possibly contam- 
inated by the Ophiuchus complex and the Orion complex, 
as suggested by Cayon & Smoot (1995). 

However, as can be seen statistically from the final three 
lines of Table 3, removal of areas possibly contaminated 
by the galactic centre and by the Ophiuchus and Orion 
complexes removes more well-matched parts of circles than 
ill- matched parts, with a increasing from 2.1 to 2.4 (the 
Clo = l,Ao = case is the one illustrated). For example, 
the central segments of panels for (i,j) — (4,6) and (5,5) 
appear consistent within about la, but both lie within 60° 
of the galactic centre. 

The positions of the Ophiuchus and Orion contamina- 
tions used are those listed as the second, third and seventh 
lines of Table 2B of Cayon & Smoot (1995). The radii of 
exclusion adopted are 57°, 9° and 29° respectively (based 
on the pixel numbers of the same table) . 

Examples of matching segments within 57° of the first 
Ophiuchus component (and not within 60° of the galactic 
centre) are the right-most segments in panels (—3, 3) and 
( — 1,6). Removal of these possible contaminating segments 
removes some conflicting segments, but also some consistent 
segments. 

Finally, if one considers the panels which are least 
affected by the galactic cut, in particular panels (—2,2), 
( — 1,1), (5,6), (6,5) and (6,6), these seem to show little 
disagreement, apart from the right-hand end of the latter 
panels. Of these, the panel ( — 1, 1) is for a pair of circles 
only separated by A = 1327h Mpc, so the agreement can 
be explained by smoothing rather than by topology. 

The right-hand end (5,6), (6,5) and (6,6) can be re- 
moved to improve the agreement by arguing that this is 
contamination by the Ophiuchus complex, as before. How- 
ever, this removal, and also removal of the galactic centre 
(as before) remove parts of the curves which do agree. If 
one takes the remaining parts of the curves, then there is 
still a broad agreement, for example if the data is further 
smoothed on a larger scale. 

But in that case, there is the possibility that what 
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Figure 2. Temperature fluctuations in COBE DMR data (the 'combined' four year Analysed Science Data Set, smoothed by 20°) around 
a selection of hypothetically identified circles in the CMB, modelled in the covering space, for f2o = 1, Ao = 0, shown against the distance 
around each circle. Thick lines show the values &T/T and thin lines show these values plus or minus the la uncertainties. The horizontal 
length is the circle circumference. (By definition, no circles arc great circles.) Each pair of circles is defined by the plane halfway between 
the observer and a topological image of the observer at iL x e x + jL y e y . The indexes i and j label each plot. The solid curves are for a 
circle at (i,j) and the dashed curves for its matched circle at (— i, —j). For the torus family of flat 3-manifolds, as in this case, the circles 
of a pair are parallel. The distance between the two members of a pair is therefore defined and is indicated here by A. Galactic latitudes 
between —20° and 20° are excluded from the plots. 



© 0000 RAS, MNRAS 000, 000-000 



COBE and Global Topology: Identified Circles 9 





Figure 3. COBE DMR 10° smoothed 'combined' ASDS map in Aitoff projection, with galactic latitude b 11 increasing from —90° to 90° 
from bottom to top and galactic longitude I 11 decreasing from 180° at the left to —180° at the right. A galactic cut of ±20° is applied. 
The positions of the three hypothesised fundamental axes determined by the three cluster images are indicated by black 'spots'. From 
right to left, these are at (—179°, 41°), determined by the vector from Coma to CL 09104+4109, (—35°, 46°), the vector from Coma to 
RX J1347. 5-1145, and at (75°, 18°) the orthogonal third axis. The positions of the clusters (Sun centred observer) are shown by white 
'spots': Coma is near the north galactic pole, the other two clusters are ~ 3° north of the corresponding axis directions. Identified circles 
for (i,j) = (—2,2) and (i,j) = (5,6) arc shown as light and dark grey curves respectively (cf Fig. 2). 



matches is merely a quadrupole or dipole component which 
has been imperfectly subtracted. The dipole subtraction for 
the ASDS's is quoted to be good but not perfect, so not only 
would a different statistic need to be calculated in order to 
measure the matching in such large scale gradients, but the 
dipole and quadrupole subtraction would have to be done 
precisely enough that this would be meaningful. 



4 DISCUSSION 

4.1 The S statistic of Cornish, Spergel & 
Starkman (1998b) 

If the error estimates are correct, then the statistics above 
are sufficient to rule out the present T 2 hypothesis. 

How successful is the correlation statistic S of Cornish, 
Spergel & Starkman (1998b) in refuting the candidate? 

The values of S shown in Table 3 are very low. Com- 
pared with Fig. 2 of Cornish, Spergel & Starkman, these 
values are quite surprisingly low. The details of the model 
for the figure are likely to differ somewhat from the present 
model, and the difference in using simulations versus obser- 
vational data (in the present case) could explain some of the 
difference. 

A value close to unity would be expected in an ideal 
case of matching circle pairs, and a value of zero for non- 
matching pairs, if the two-point auto-correlation function of 
the temperature fluctuations were zero. In reality, there is 
an intrinsic non-zero auto-correlation function, so a value 
showing this weak a correlation is unexpected. 

If standard arguments implying the relative weakness of 




100 200 300 

mod( l" - 74.5°, 360° ) 

Figure 4. Dependence of cr, i.e. the standard deviation from iden- 
tity of temperature fluctuations around hypothetically multiply 
imaged circles on the CMB, in units of the observational uncer- 
tainties [eq. (1)]. The value at galactic longitude I 11 = 74.5° is 
for the Roukema & Edge (1997) T 2 cluster-based hypothesis and 
is clearly greater than a = 1. Points at other longitudes are for 
other hypotheses equivalent to that of Roukema & Edge, except 
for a rotation in galactic longitude. The latter are unmotivated 
by observations. This and the following three figures are for the 
20° smoothed 'combined' map. 



the Doppler effect are invalid, could this be due to a domi- 
nant contribution by the Doppler effect, which would cause 
genuinely anti-correlated temperature values to cancel when 
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Figure 5. The mean temperature difference d [eq. (2)] along iden- 
tified circles is plotted, as for Fig. 4. 




100 200 300 

mod( l 11 - 74.5°, 360° ) 



Figure 6. Dependence of S [eq. (3)], the correlation statistic 
defined by Cornish, Spcrgcl & Starkman (1998b), along identified 
circles as for Fig. 4. 



evaluating the 5* statistic? The parameter 5*£) pp does not 
support this. Its value is smaller in absolute value than that 
of S, presumably because typical values of fii-fij are smaller 
than unity and there is no underlying ant i- correlation. 

Figures 4-7 show that the small value of S is just a co- 
incidence. These plots show equivalent hypotheses to that 
based on the three clusters, apart from a rotation in galac- 
tic longitude. This avoids any systematic effects due to the 
galactic cut or to other functions of galactic latitude. 

Assuming that the hypothesis is false, the ranges of val- 
ues of the statistics for different longitudes show typical val- 
ues of the parameters for the T 2 models equivalent to the 
hypothesised one (equivalent apart from the lack of obser- 
vational motivation in the rotated models). 

Fig. 6 shows why the low value of S for the Roukema 
& Edge (1997) candidate is just a coincidence. The position 
of the candidate just happens to lie at a longitude where 
both d and S have low (absolute) values. The longitudi- 




100 200 300 

mod( l 11 - 74.5°, 360° ) 

Figure 7. Dependence of »Srj pp [ ec l- (4)]j a modification of the S 
statistic suggested for detecting correlations of the Dopplcr effect, 
along identified circles as for Fig. 4. 



nal dependence of d clearly has a sinusoidal-like component, 
which could be explained by some large scale feature in the 
temperature map. The longitudinal dependence of S shows 
more complexity, and provides a distribution of maximum 
values about a factor of two to four times lower than those 
of '5 m ax' in Cornish, Spergel & Starkman's simulation. 



4.2 The 5*£) pp statistic 

The shape of the Sdopp distribution (Fig. 7) is quite differ- 
ent from that of the S distribution (Fig. 6). Typical values 
of SrjQpp are roughly half that of the typical S values. An 
order of magnitude estimate to explain this is that (very) 
roughly half of the pairs of matched pixel vectors are roughly 
orthogonal, and so contribute very little to SiWp. 

Some features in the S distribution, e.g. several local 
maxima and the global minimum in Fig. 6, can be identified 
with similar features in Fig. 7. These can be interpreted as 
features due to large matched circles for which the matched 
pixel vectors are roughly parallel. 

Apart from the 'missing' contribution of orthogonal vec- 
tors in Fig. 7, the differences between the two figures also 
come from roughly opposite vector pairs, i.e. from pairs of 
small circles for large i and j values, where the 'second ob- 
server' is close to 2rgLg from the first observer. 

This demonstrates an approach which will be partic- 
ularly useful for the study of MAP and Planck Surveyor 
data. The statistic S-£) p p as denned in eq. (4) will provide 
a simple method of searching for roughly opposite, small 
identified circles in which a Doppler component is strongly 
expected to make a significant contribution. 

If a candidate 3-manifold is detected in a new CMB 
map, then comparison of 'false' hypotheses rotated in galac- 
tic longitude such as here will provide a simple way of check- 
ing the significance of the detection, avoiding the need for 
simulations. 
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4.3 Comparison of COBE data to arbitrary T 2 

models via simulated temperature fluctuations 

In the absence of specific candidates for torus models, several 
authors have previously compared COBE data to sub-classes 
of these models, but using the perturbation simulational ap- 
proach which requires assumptions on perturbation statis- 
tics as described above (Sect. 1.2). The closest work among 
these to the present study is that of de Oliveira-Costa et al. 
(1996), who considered T 1 and T 2 models. 

Because of the presence of one (or two) long axis(es) , the 
'large-scale cutoff in Fourier mode' argument for T 3 models, 
which can be valid for the assumptions of Gaussian random 
amplitudes in the Fourier modes, does not apply in these 
cases. 

These authors defined a statistic, So, which is similar 
to a, but which includes multiply imaged pixel pairs along 
with many non-multiply imaged pixel pairs. The advantage 
of this is that the statistic is a minimum value, for which the 
orientation of the fundamental polyhedron is not assumed. 
The disadvantage is that non-multiply imaged pixel pairs 
are indiscriminately mixed with multiply imaged pixel pairs. 
The points (a), (b), (d), (e) and (f) (Sect. 1.2.1) also apply 
to this study. 

This enabled simulations, based on assumptions regard- 
ing the power spectrum of density perturbations, to be made 
in a parameter space comprising only the lengths of the 
fundamental polyhedron and two random variables for the 
Fourier modes, but not the orientation of the fundamen- 
tal polyhedron. The statistic could therefore be calculated 
within practical computing capabilities. 

Although this So parameter is not identical to a, and 
it includes noise due to non-corresponding pixels, the value 
obtained is 2.59, not too much larger than the value of a 
which only compares matching pixels. This suggests that the 
So parameter could still be useful for application to Planck 
and MAP data. 



5 CONCLUSIONS 

The COBE Advanced Science Data Sets (ASDS) have been 
used to try to refute the candidate for the global topology 
of the Universe defined by supposing that the three clusters 
Coma, RX J1347.5-1145 and CL 09104+4109 are three topo- 
logical images of a single cluster and determine two genera- 
tors of the 3-manifold, and that the third generator is large 
and orthogonal to the first two. The candidate is a T 2 x R 
candidate, for a flat, i.e. Qo + Ao = 1 universe, where Ao = 0. 
Given the assumptions that 

(i) the error estimates given in the ASDS's are sufficient as 
good estimates of the total random plus systematic errors, 
and 

(ii) the temperature fluctuations measured are dominated 
by the naive Sachs- Wolfe effect which averaged over spheri- 
cal sections results in locally isotropic emission, 

then, 

(i) using either the 'subtracted' or the 'combined' ASDS, 

(ii) with smoothing at either 10° or 20°, and 

(iii) with or without making cuts for possible contamina- 



tion by the galactic centre and by the Ophiuchus and Orion 
complexes, 

the difference statistics a and d clearly refute the T 2 hy- 
pothesis to more than 94% significance. (This is estimated 
by a Kolmogorov-Smirnov test comparing the list of tem- 
perature differences to a Gaussian distribution of mean zero 
and standard deviation unity.) 

A galactic cut of ±20° is used for all the above options. 

The subtracted set more strongly rejects the hypoth- 
esis than the combined set. Smoothing at 20° leads to a 
stronger rejection than for smoothing at 10° . Removal of the 
suspected contaminating regions strengthens the rejection, 
suggesting that these introduce more spurious well-matched 
portions of circles than ill-matched portions. 

This result does not rely on assumptions on the form 
of the fluctuation spectrum or on the distributions of ampli- 
tudes and phases of the spectrum. 

The validity of assumption (i) will be tested by the in- 
creasingly numerous small angle measurements and by the 
Planck and MAP missions. 

In the case of a non-zero cosmological constant (Ao = 0, 
e.g. Fukugita et al. 1990), the T 2 candidate has not yet been 
ruled out, since a modelling of the integrated Sachs- Wolfe 
(ISW) effect would be required. In this case, assumption 
(ii) would be invalid to the extent that the ISW contributes 
to the COBE CMB (Sect. 7.3, White, Scott & Silk 1994). If 
the magnitude of this component is of order of the estimated 
noise, then this would be equivalent to assumption (i) failing 
due to systematic error. 

It is interesting to note that the matched circles prin- 
ciple provides a relatively clean and optimal method of re- 
futing specific candidates for global topology. For data of 
the COBE quality, it is clear that there is strong sensitivity 
to the validity of the noise estimates. Since this is the case 
for the matched circles method, it must be an even stronger 
caveat for the perturbation simulation methods, which have 
mostly been applied to T 3 models for which L x — L y — L z , 
for which the generators have been assumed to be mutu- 
ally orthogonal, and for which simulations of the fluctuation 
spectra were required. 

Since the effectiveness of this new method of applying 
the identified circles principle in order to directly refute a 3- 
manifold candidate (apart from the caveats) has been shown, 
future work would be to use the same method to refute cer- 
tain classes of flat (or hyperbolic) multiply connected uni- 
verse 3-manifolds with COBE data. This would complement 
previous work, which relied on assumptions on the pertur- 
bation statistics which are probably inconsistent with the 
hypotheses which were tested. 
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